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Abstract. We demonstrate the benefits of a convex geometric perspective for 
questions on chemical stoichiometry. We show that the balancing of chemical 
equations, the use of "mixtures" to explain multiple stoichiometry, and the 
half-reaction for balancing redox actions all yield nice convex geometric inter- 
pretations. We also relate some natural questions on reaction mechanisms with 
the enumeration of lattice points in polytopes. Lastly, it is known that a given 
reaction mechanism imposes linear constraints on observed stoichiometries. 
We consider the inverse question of deducing reaction mechanism consistent 
with a given set of linear stoichiometric restrictions. 
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A common question encountered in chemistry is the balancing of a chemical 
equation, and it has been widely discussed in the chemical education literature (see 
for example the review article by Herndon|12|). It is well-known in the literature 
(though perhaps not among chemistry students) that the question admits a linear 
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algebraic formulation. Here we examine an aspect previously unexplored, namely 
a convex geometric approach. We show that after specifying the species (i.e. the 
atoms or compounds involved in a reaction) on each side of a chemical equation, 
convex polytopes (a bounded intersection of half-spaces) may be used to provide 
a visual illustration for questions of existence and uniqueness of balancings, and 
though generally applicable, it is particularly effective for chemical equations in- 
volving only neutral species where no more than three or four elements in total are 
represented. More generally, we explore the utility of convex geometric representa- 
tions in chemical stoichiomctry, and this forms the principal theme of the current 
article. 

The convex geometric approach complements well the usual linear algebraic 
method which we summarize in Section[TJ with particular attention to issues regard- 
ing existence and uniqueness of solutions. That section concludes with a geometric 
interpretation for the result of the computation. Section [2] discusses the aforemen- 
tioned geometric approach to the balancing question and its connection with the 
algebraic approach. The section includes various examples with diagrams illustrat- 
ing the potential visual appeal and pedagogical value of this approach. In Section [3] 
we continue to emphasize the utility of the geometric approach by explaining vari- 
ous chemical practices in this geometric framework. In particular we show that the 
practice of explaining chemical equations arising in contexts of multiple balances as 
mixtures of equations with unique balances is mathematically justified. Section [4] 
discusses adaptations of the geometric approach in the presence of charged species. 
There is an extended discussion of the "half-reaction method" commonly used to 
balance oxidation-reduction reactions. The geometric approach nicely elucidates 
the scope and limits of this practice. Section[5]turns attention from individual equa- 
tions to reaction mechanisms. There are two major considerations in this section: 
First, we show how natural existence and enumerative questions about mechanisms 
lead to the much-studied problem of counting lattice points in polytopes. Second, 
we investigate the moduli of reaction mechanisms consistent with a given overall 
reaction and observed stoichiometric linear dependencies. The latter complements 
the "forward" analysis of Missen and Smith [TS] in deducing stoichiometric linear 
dependencies given a reaction mechanism for an overall reaction. Finally, the Ap- 
pendix (Section [6]) includes a geometric approach to some of the observations in 
Section |1 

Below are some highlights of the results and perspectives of this paper: 



(Section [2]) Questions of existence and uniqueness of balancings for a chem- 
ical equation may be phrased in terms of the intersections of convex poly- 
topes. This convex geometric approach is visually appealing, perhaps peda- 
gogically beneficial, and intimately connected to the known linear algebraic 
formulation of the problem. 

(Proposition |3.7[ ) If a chemical equation admits multiple balancings, then 
any such balance may be realized as a "mixture" of reactions admitting 
unique balances. 



(Proposition 4.111 If a redox reaction admits a balance, then the "half- 



reaction method" will generate a balance for the equation. However, in 
case of multiple balances, not all possible balances may be obtainable via 
the "half-reaction" method. 
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• (Remark |5.3| ) In cases of multiple balances, minimizing the sum of coeffi- 
cients in the balance may not uniquely determine a balance. 

• (Subsection 5.3) It is known that any given mechanism imposes linear de- 
pendencies on observed quantities of species. The "inverse" problem of 
deducing mechanisms consistent with given observed linear dependencies 
on species yields a collection of mechanisms parametrized by subspaces 
that can be explicitly identified given knowledge of reaction intermediates. 
This does not completely solve the inverse problem since convex geometric 
and order considerations still need to be imposed. The latter is resolved 
but the former may be computationally expensive. 

Computations on convex polytopes were done with Matthias Franz's Maple pack- 
age Convex [TU] with the exception of lattice point enumeration which was com- 
puted using Verdoolaege's program barvinok 24J. 



1. The Algebraic Approach 

The first subsection summarizes the linear algebraic approach, while the second 
geometrically interprets the result of the computation. 

1.1. The Algebraic Approach in a Nutshell. We begin with an example il- 
lustrating the connection between balancing chemical equations and linear algebra. 
Further references and examples are provided by Missen and Smith IB] , Blakley^Q 
and AlbertypQ among many others. 

1.1. Example. Suppose we want to know all possible ways to balance a reaction 
involving hypothetical neutral species XY, YZ, XYZ2 where X, Y, Z are distinct 
elements. The algebraic method associates to each species a 3-vector, in this case 




vyz = y. j VXYZv 

where we have tacitly ordered the elements X, Y, Z so that the first component 
indicates the number of atoms of X, the second the number of Y atoms, etc.. 
Mass-conservation is then reflected by finding rational numbers ai such that 

dlVXY + CL2VYZ + a-iVxYZ 2 = 

One just multiplies by the least common denominator of the fractions to obtain 
an integral solution. Thus, if a balancing exists, then a matrix M with the above 
vectors as its columns has a non-trivial nullspace NS(M): that is, the equation 
Mx = does not have just the vector = (0, . . . , 0) T as a solution. For example, 
we may set 

A 

M= 1 

\o 

(The ordering of the columns does not affect non-triviality of the nullspace.) One 
computes that NS(M) = so that no balancing exists. Conversely, if NS{M) 7^ 0, 
then since M has integer entries, one may find a basis for NS(M) having only 
rational entries. Hence, a balancing exists iff NS(M) ^ and it is furthermore 
unique (up to multiplicative factor) precisely if the nullspace is one-dimensional. 
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We formalize the above as follows: Suppose we have m neutral species involving 
a total of n elements which we place in some order. To each species we associate a 
n-vector where the k-th component is the multiplicity of the fc-th element in that 
species. Denoting these species vectors by Vj, we seek rational numbers <2j such 
that 



(1.2) J2 a ^ = ° 



»=i 

Or equivalent, let M = (vj . . . v m ) beanxm matrix with the vectors Vj as columns. 
Then, a balancing exists iff dim(NS(M)) > and the balancing is unique (up to 
multiplicative factor) iff dim(NS(M)) = 1. 

Note that the method does not indicate which species are reactants and which 
products. Only mass-conservation is being enforced and no directionality is implied, 



though for a given solution to Equation (1.2 1 the species are separated into two 
collections depending on the sign of the coefficients a* . (Those for which <Zj = are 
not involved in the balancing.) For convenience, we will refer to these two groups 
as "reactants" and "products" without implying any directionality. 

Now, suppose we stipulate that certain species are "reactants" and the remaining 
"products." In the case where dim(NS(M)) = 1, the balancing is unique and the 
species are partitioned into two groups uniquely. However, for dim(NS(M)) > 
1, certain partitions may not be realizable and even if realizable, the balancing 
may not be unique (even up to multiplicative factor). Mathematically, we seek 



rational coefficients at in Equation (1.2) such that the signs of those associated to 
"reactants" are opposite of those associated to "products." Computationally, one 
may proceed as illustrated in Blakley[4] by computing a basis {b,} for the nullspace. 
Since now the components of bj index the species, we seek linear combinations of 
hj such that components indexing "reactants" have signs opposite those indexing 
"products." 

1.3. Example. Suppose we have species XY, XZ, YZ, XYZ, X5Y5Z2 partitioned 
into "reactants" {X,Y,XYZ} and "products" {XZ,YZ, Z 5 Y 5 Z 2 }. Ordering the 
elements X, Y, Z in that order, and defining 

/l 1 1 5> 
(1.4) M = 1 1 1 5 

\0 1 1 1 2y 

we compute that dim(NS(M)) = 3 with the vectors 
(1.5) 

bi = (0,1, -1,1,0, Of b 2 = (1,0, -1,0,1, Of b 3 = (-3,-3,-2,0 ! 0,lf 

as one possible basis for the nullspace. Note that the components of b^ index the 
species X, Y, XYZ, XZ, YZ, Z 5 Y 5 Z 2 in this order. Choosing the "reactants" to have 
non-positive coefficients, we thus seek linear combinations cibi + C2b2 + C3b3 such 
that the first three components are non-positive and the remaining non-negative: 

c 2 - 3c 3 < ci > 

ci - 3c 3 < c 2 > 

-ci - c 2 - 2c 3 < c 3 > 
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We further require that all resulting six components are rational numbers (see fol- 
lowing discussion) . A geometric approach for this balancing is presented in Example 
2.51 below. 



We may formalize the preceding as follows: Suppose we partition the m neutral 
species into r "reactants" and p "products." After possible re-ordering, we may 
assume vectors vi, . . . , v r correspond to the "reactants" and vy+i, . . . , v m to the 
"products," where r+p = m. Without loss of generality we may choose "reactants" 
to have non-positive coefficients. Then, each balancing corresponds to a collection 
of rational numbers a, such that $^ i=1 ctiV.j = now with the added stipulation 
that ai, . . . , a r be non-positive and a r +i, . . . , a m non-negative. Let M — (vi . . . v m ) 
and {hj} be a basis for the nullspace. Since M is a matrix of integers, we may 
assume that the components of the basis vectors hj are all rational. Then, the 
desired rational Oj correspond precisely to rational linear combinations since the 
hj constitute a basis for NS(M) and hence are linearly independent: 
(1.6) 

' a i\ k f c i\ 

= c j b j = ( b i • • • b fc) : ai,...,a r <0 a r+x , . . . , a m > 

w 

where k = dim(NS(M)) and the Cj are rational. 

1.2. A Geometric Interpretation for the Result of the Algebraic Method. 

We provide a geometric interpretation for the system of inequalities in Equation 



(1.6). Each component ai defines a linear inequality in the variables Cj. Geomet- 
rically, this represents a half-space in K fc passing through the origin, where recall 
k = dim(NS(M)). Thus, simultaneous solutions to the system of inequalities in 
Equation (1.6) are geometrically represented by the intersection locus Q of m half- 
spaces in R*. Such an object is called a polyhedron in K fc . Since we want rational Cj, 
the possibilities are precisely the rational points (i.e. points where all coordinates 
are rational numbers) contained in the polyhedron where two such rational points 
represent the same balancing (up to multiplicative factor) if they are on the same 
line through the origin. Thus, the "moduli" of balancings with given "reactants" 
and "products" is realizable as the image of Q C ]R fc under its projectivization, 
hence a subset of rational projective space Pq. We will return to a discussion of Q 
near the end of the next section. In summary, we have a geometric representation 
for all possible balancings of a chemical equation when species for both sides are 
initially specified. 

1.7. Example. Consider the oxidation of nitric oxide (NO) to nitrogen dioxide 
(A r 2 )- The reactants are NO,0^ and the products are N02,C>2- Note that 
this example involves two allotropes of oxygen. We order the elements by O, N in 
that order. Defining 

'1 3 2 2^ 



( L8 ) M =n o i o 

we compute that dim(NS(M)) = 2 with the vectors 

(1.9) bx = (0,-2,0, 3) T b 2 = (-3,-l,3,0) T 

as one possible basis for the nullspace. Note that the components of b^ index the 
species NO, 3 ,N0 2 , 2 in this order. Choosing the reactants to have non-positive 
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coefficients, we thus seek linear combinations cibi + c 2 b2 such that the first two 
components are non-positive and the remaining non-nonegative: 

-3c 2 < 3c 2 > 

-2c x - c 2 < 3ci > 

The resulting polyhedron Q is the set {{c\,C2)\ci > 0} and each rational point 
within Q indexes a balancing. A geometric approach for this balancing is presented 
in Example |2.7| below. 

2. Geometric Approach to Balancing 

Now we examine the balancing question anew but guided by geometric con- 
siderations. Suppose as before we partition m neutral species into r "reactants" 
and p "products" and that these are represented by the collections Vi , . . . , v r and 
v r+ i, . . . , v m respectively, where r + p = m. Then, balancings are precisely equali- 
ties of a non-negative rational linear combination of the former set of vectors with 
a non-negative rational linear combination of the latter set: 

r m 

(2.1) Y. riVi = H Pj' v i 

i—l j—r+1 

where the coefficients r^Pj > are rational. 

Recall that the set of all non-negative scalings of a vector is geometrically rep- 
resented by a ray based at the origin in the direction of the vector. Non-negative 
linear combinations of a set of vectors determine a mathematical object called the 
(polyhedral) cone spanned by the vectors with the origin as the vertex of the cone. 
Since we only care about the rational numbers ri,pj up to multiplicative factor, we 
slice both cones simultaneously by a hyperplane intersecting all the positive x j-axes 
for i = 1, ... This hyperplane then intersects each cone transversely and the 
intersection loci are bounded polyhedra, i.e. (convex) polytopes, which we may call 
the "reactant" and "product" polytopes. Indeed, each is the convex hull of the points 
given by the intersection of the various defining rays (either those associated with 
"reactants" {v 1; . . . , v r } or those with "products" {v r+1 , . . . , v TO }) with the slicing 
hyperplane. Without loss of generality, we may assume the hyperplane is described 
by an equation of the form n ■ x = h where n, h consist entirely of rational numbers. 
Then, given a choice of "reactants" and "products," a necessary condition for the 
existence of a balancing is that the intersection of these two polytopes is non-empty. 
The intersection is another polytope, which we will call the intersection polytope, 
and if non-empty, it necessarily contains a rational point. Conversely, one can show 
that each rational point in the intersection polytope corresponds to some balancing 
of a chemical equation with the specified grouping of species. (See the discussion 
at the end of this section or for a geometric argument see the appendix.) 

Since we assume that the coefficients ri,pj > are non- negative, all the compo- 
nents of the scaled vectors riVi,pjVj are non-negative, and a canonical choice for 
the hyperplane is X)T=i Xl = Notice that this hyperplane intersects each positive 
x^-axis for i = 1, . ..,n as required. The intersection of this canonical hyperplane 
with the region {(xx, . . . ,x n )\xi > 0} is a (n — l)-dimensional "triangle," called a 
(n— 1)- simplex. The coordinates (x±, . . . , x n ) sum to unity and are called barycentric 
coordinates. Each species vector Vj corresponds to a point with barycentric coordi- 
nates specifying the proportions of each element within the species. For example, if 
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Figure 1. Case where no balancing exists for any grouping of 
reactants and products 



X, Y, Z are the only elements involved and ordered thus, then the barycentric coor- 
dinates for the species X2Y3Z4 is (2/9, 1/3,4/9). Taking a non-chemical example, 
the RGB-system of colors describes colors in terms of proportions (r, g, b) of red, 
green, and blue, respectively, where r + g + b = 1. 

Note that rational points on the boundary of the intersection polytope corre- 
spond to balancings in which not every species is present. Hence, if we require that 
all species be present, then if an equation can be balanced then it can be balanced 
in infinitely many distinct ways unless the intersection polytope is a single point. 
However, this is only a necessary and not a sufficient condition for the uniqueness 
of balance. We will see later that additionally we need the generating vectors for 



the reactant and product cones each to be linearly independent (see Equation 2.18 1 . 
For example, this is satisfied if the reactant and product polytopes are line segments 
meeting at a point in the relative interior of each. 

We now illustrate this perspective with several examples: 

2.2. Example (Non- Existence of Balancing) . Suppose we have "reactants" {XY,YZ} 
and "product" {XYZ 2 } where X, Y,Z are distinct elements and all species are 
neutral. If we order the elements as X, Y, Z, then the barycentric coordinates of 
the "reactant" species are (1/2, 1/2, 0) T and (0, 1/2, 1/2) T , and that of the "prod- 
uct" species is (1/4, 1/4, 1/2) T . The "reactant" polytope is the convex hull of the 
barycentric coordinates for the "reactant species," namely a line segment joining 
(1/2,1/2, 0) T and (0,1/2, 1/2) T whereas the "product polytope" is just a single 
point. See Figure [T] Since the line segment and point are disjoint, we conclude 
that no balancing exists for that particular choice of "reactants" and "products." 
In fact, from the diagram we can easily see that there is no way to partition the 
species into "reactants" and "products" to obtain a situation where a balancing 
exists. This concurs with the conclusion in Example 1 1 . 1| obtained via the algebraic 
method. 
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2.3. Example (Unique and Non-Unique Balancings). Consider the neutral chemical 
species H 2 ,H 2 0, CH i: C0 2 , CO. Ordering the elements H, O, C in that order, the 
barycentric coordinates are respectively, 

(2.4) 

See Figure [2] As evident from the diagram, to obtain a non-empty intersection 
polytope, we need to have at least two "reactants" and two "products" each. By 
the comments above, we have an unique balancing precisely when four of the five 
species are involved and are grouped pairwise, and the lines so determined intersect 
uniquely One easily checks from the diagram that this yields the groupings: 










{CH 4 ,C0 2 } 
{H 2 ,CO} 
{H 2 ,CO} 
{H 2 0,CO} 
{H 2 0,CO} 


{H 2 ,CO} 
{CH Al C0 2 } 
{CH4,H 2 0} 
{CH Al C0 2 } 

{H 2 ,C0 2 } 



(Compare with cases (a)-(e) in Example 6 of Blakley[3]). 



2.5. Example (Non- Unique Balancing). Consider "reactants" {X,Y,XYZ} and 
"products" {XZ,YZ,X 5 Y 5 Z 2 } where X,Y,Z are distinct elements ordered thus 
as in Example |1.3| The barycentric coordinates are 
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Figure 3. Case where infinitely-many distinct balancings exist for 
given choice of reactants and products 



O, , 0, NO 



(1,0) NOj (0,1) 

Figure 4. Case where allotropes are involved 

respectively. See Figure [3] The intersection polytope is a quadrilateral, and hence 
there are infinitely many distinct balancings under this grouping of species. For 
example, the rational points (3/8, 3/8, 1/4) T and (2/5, 2/5, 1/5) T are in the interior 
of the intersection polytope with 

1 





1 

5 

yielding distinct balancings: 

2X + 2Y + AXYZ = XZ + YZ + X 5 Y 5 Z 2 
8X + 8Y + 8XYZ = XZ + YZ + 3X 5 Y 5 Z 2 

respectively. 



2.7. Example (Allotropes). Consider the oxidation of nitric oxide (NO) to nitrogen 
dioxide (N0 2 ) as discussed earlier in Example 1.7 The reactants are NO, O3 and 
the products are N0 2 ,0 2 . Ordering the elements by 0,N in that order, both 
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allotropes are represented by the barycentric coordinates (1,0) T whereas those for 
NO,N02 are (1/2,1/2) T and (2/3, 1/3) T respectively. The intersection polytope 
is a line segment {(1 - t,t) T \0 < t < 1/3}. See Figure [I| Writing (1 - t,t) T for a 
point in the intersection polytope, we have the representations: 

w> (v)-K$) + *-»>®-'(D +i ^G) 

M ( 1 7 , )-«P +(1 - i "(3-'(0 + L r ! (?) 

(2.10) 

relative the react ant and product cones, respectively. Thus, we have a one-parameter 
family of distinct balancings: 

1 — 2/ 1—3/ 
(2.11) /iVO + — ^— 3 = /7V0 2 + — ^— 2 

indexed by rational t such that < / < 1/3. If we insist that all reactants and 
products be present, then we have strict inequalities. For example, setting t = 1/4 
and clearing denominators yields the balancing: 6N O + 4C*3 = 6NO2 + 3C>2 • The 
case where t = 1/5 corresponds to the usual equation for the oxidation: NO + O3 — > 
N0 2 + 2 . 

2.12. Remark. In practice, one knows the relative amounts of reactants and prod- 
ucts. Exact knowledge of either determines an unique ray in the corresponding 
cone. Hence, if a solution exists, this information is sufficient to identify an unique 
point in the intersection polytope. Because of measurement error and uncertainty, 
we have instead a neighborhood of a ray (mathematically, a neighborhood of a 
point in projective space), and hence possibly infinitely-many balancings. In prac- 
tice, chemists choose balancings with "small" coefficients, though a mathematical 
formulation for this "rule-of-thumb" is unclear. See also Remark 15.31 . 

2.13. Remark. Note that we can scale the intersection polytope so that it has 
integer vertices. Let d be its dimension and write v(n, d) for the number of points 
in the intersection polytope with denominator (in lowest terms) no greater than n. 
Analogously, let i>o(n, d) denote the number of such in the interior of the polytope. 
By convention, we set v(0,d) = 1 = ^o(0, d). It is known that both v,vq are 
polynomials of degree d in n (see Stanley |21j). In particular, knowledge of d 
other values determine either v or vq completely. For example, we can scale the 
quadrilateral that is the intersection polytope in Figure [3] to have integer vertices 

(15, 15,6) T (16, 10, 10) T (10, 16, 10) T (12, 12, 12) T 

Then, with regard to this polytope, we can count that in its interior there are 
16 integer points and 33 points with denominator at most two so that vq(ti,2) = 
n 2 + 14n + 1. However, relating denominators to the least integral coefficients of 
a balancing is not direct, as Example |2.7| shows. We will return to enumerative 
questions in Subsection |5.2| 

Now, we describe the connection between the geometric and algebraic approaches: 
Recall that the algebraic approach parametrizes possible balancings by rational 
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points of an unbounded polyhedron Q in R k . Consider the linear transformations 

(2.14) B: Q (b l^ fc) r 

(2.15) C r :R r + {vi -^> r) R n 

(2.16) C p :R p + iVr+ ^ Vm) R n 

where R*j_ denote non-negative values and we have interpreted matrices as linear 
maps. Define R = C r o ir r and P = C p o (— 7r p ) where 7r r (resp. 7r p ) denotes 
projection onto the first r (resp. last p) coordinates in R m . Then, B : Q — » M m 
is the centralizer of the maps i?, P : R m — > R n (in particular, they coincide on the 
image of B), and RB, PB : Q — > R n each map the polyhedron Q in a linear fashion 
onto the intersection cone. These maps fit into a commutative diagram: 

(2.17) 




where R = Tt r o B and P — ~ir p o B. 

Let 1Z, V be the reactant and product cones, respectively and X = TZ C\ V the 
intersection cone. Note that 1Z = imC r and V = imC p . Define the matrix M — 
(vi . . . v m ) whose columns are those of C r ,C p respectively. Note that NS(M) = 
im(B) so that dim(i\T5(M)) = dim(Q) the dimension of the unbounded polyhedron 
in the algebraic approach. By a fact from linear algebra (see Section [6]), one has 

(2.18) dim(Q) = dim(NS(M)) = dim (I) + dim(ker(CV)) + dim(ker(C p )) 

so that of the polyhedron Q from the algebraic approach is not in general linearly 
isomorphic to the intersection cone. It is so iff C r , Cp are injective, that is, the reac- 
tant and product species vectors are each linearly independent sets. One example 
where this fails is a reaction involving either more reactant or product species than 
elements involved. 

Intuitively, the nullspace NS(M) measures linear relations among the reactant 
and product species vectors. These may be among reactant vectors only, among 
products vectors only, or among vectors drawn from both reactant and product 
species, hence the three terms in Equation |2.18| Letting s,r,p be the number of 
species, reactants, and products respectively, we have in species space S = R s 
the (s — l)-simplex A s_1 = A r_1 * A p_1 where A 1 "" 1 , A^" 1 are mutual cofaces 
defined by the reactant and product vectors, respectively. The species vectors 
define partial maps A r ~ 1 , A p_1 =^ £ whose images are the reactant and product 
polytopes, respectively. The bijectively image of the moduli polyhedron Q in species 
space S is possibly "partially collapsed" under the "elemental proportion" map M 
into the intersection cone in element space £ due to possible linear dependencies 
among generating rays for the reactant cone and also for the product cone. In 
summary, under the geometric approach, we gain often a visual determination of 
existence issues at the expense of possibly losing count of the algebraic dimension 
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of possible multiple balancings and hence a slightly more involved criterion for 
uniqueness. 

3. Polytopal Explanation for Some Chemical Practices 

In this section we examine several common practices in chemical stoichiometry 
and show that they all have simple polytopal explanations. 

3.1. Inspection Methods. In the chemical education literature, there are numer- 
ous "inspection methods" for balancing chemical equations (for example [23|.[15|). 
Mathematically, these often amount to establishing an ad hoc elimination order 
on the elements, or in geometric terms, matching the projections of the reactant 
and product polytopes onto various axial directions. For example, the first two 
of "Ling's Rules for Balancing Redox Equations by Inspection" as presented in 
Kolb[H] are: 

• "Step 1: Locate any elements that must have the same coefficient in the 
balanced equation, those appearing only once on each side of the equation 
and in equal numbers on both sides. Mark these terms with arrows." 

• "Step 2: Locate any elements that appear only once on each side of the 
equation but have unequal numbers of atoms. Balance these elements first." 

This is illustrated in the following example: 

3.1. Example. |15j In balancing 

S + HN0 3 -> S0 2 +NO + H 2 

by step 1, we see that S,S02 have identical coefficents and so also HNOs,NO. 
By step 2, we balance hydrogen, which forces the balance for nitrogen, then sulfur, 
then oxygen yielding: 



Kolb's article [15] even presents a case where Ling's method does not work, and 
she resorts to linear algebra, the proper algebraic framework. 

3.2. Stoichiometric Restriction. Sometimes one knows from experiment that 
certain species react in particular ratios, usually from kinetic and mechanistic con- 
siderations (see Section [5] for details). Missen and Smith [17] provide an example 
involving a permanganate-peroxide reaction in acidified aqueous solution: 

H 2 SO A + KMnSOi + H 2 2 -> 2 + H 2 + K 2 SO± + MnSO A (unbalanced) 

One checks (for example via Convex[1G\ or Polymake[TT]) that the intersection 
polytope is 1-dimcnsional. However, it is known from experiment that KMnOi 
reacts with H 2 2 in a 2 : 5 ratio. Algebraically, this may be reflected by augmenting 
an additional row to the matrix M mapping species to elements. For example, 
ordering the species as they appear left-right in the above equation, the augmented 



Geometrically, we incorporate the restriction by replacing the generating vertices 
VKMnS0 4 ,v H2 o 2 with a linear combination (2/7)v KMnS o 4 + (5/7)v H2 o 2 . One can 
compute that this alteration reduces the dimension of the reactant polytope. Either 




matrix is 




M 

5 
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Figure 5. Reaction as Mixture of Reactions 

way, we compute that the intersection polytope is now O-dimensional so that there 
is an unique balance under this stoichiometric restriction. Similarly, if is the 
generating vertex indexed by the z-th specie, and it is known experimentally that 
reactants Ri 1 , . . . , Ri k react in a ratio of : • • • : , then we replace the associated 
generating vertices with the single vertex given as a weighted linear combination: 

(3-2) v = ^v l3 r = 5>«, 

3=1 3=1 

A similar observation applies if certain products are know to be produced in par- 
ticular ratios. The intersection from the resulting polytopes may or may not be of 
lower dimension. 

More generally, if there is linear relation among the species, it defines a hyper- 
plane H that interesects the simplex A s_1 C S. Intersecting the image M(H n 
A" -1 ) C £ with the intersection polytope then yields the new modified intersection 
polytope. Algebraically, one just adds these relations as additional rows to the 
matrix defining the elemental composition map. 

3.3. Mixtures of Reactions. In the chemical education literature, reactions ad- 
mitting multiple balancings are often explained as a mixture of several other re- 
actions. This is tantamount to writing a balance as a rational combination of 
balancings relative subpolytopes of the reaction and product polytopes. 

3.3. Example. [Mj Kolb explains the reaction: 

3HC10 3 HC10 4 + Cl 2 + 20 2 + H 2 
as a mixture of two reactions: 

7HC10 3 bHClOi + Cl 2 + H 2 

4HC10 3 2Cl 2 + 50 2 + 2H 2 



14 



JER-CHIN (LUKE) CHUANG 




Figure 6. Reaction as Mixture of Reactions 

where the first plus twice the second yields the desired overall reaction. From the 
polytopal diagram (see Figure[5]), we see that HCIO3 is in the interior of the product 
polytope, and hence the reaction admits multiple balancings. The first reaction 
expresses the reactant via the subpolytope spanned by products HCIO4, Cl 2 , H 2 
whereas the second via that spanned by products Cl 2 ,0 2 , H%0. 

3.4. Example. [Mj Kolb also explains the reaction: 

350 2 + 7C -> CS 2 + S + 6CO 
as a mixture of the reactions 

S0 2 + 2C -> S + 2CO 

2S0 2 + 5C -> CS 2 + ACO 

Again, from the polytopal diagram (see Figure [6]), the existence multiple balancings 
is clear. Note that as in the preceding example, the particular choices for component 
reactions each admit unique balancings. 

3.5. Example. |14j In the equation 

NaCl0 2 + NaCIO -> NaCl0 3 + NaCl 
chlorine appears in four different oxidation states. Hence, application of the usual 



half-reaction method is somewhat confusing (see discussion in Subsection 4.3 be- 
low). Thankfully, the equation as written is already balanced, though it also admits 
multiple balancings. From the polytopal diagram (see Figure [7|, we see that the 
species are actually collinear, since each has Na, CI in a 1 : 1 ratio. Kolb explains 
this example as a combination of the reactions 

3NaCl0 2 -> 2NaCl0 3 + NaCl 
SNaCIO -> NaCl0 3 + 2NaCl 
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Figure 7. Reaction as Mixture of Reactions 



In all the preceding examples, the choice of component reactions was chemically- 
motivated but also happened to admit unique balancings. Indeed, one may say they 
were chosen in part to be uniquely determined, since the paradigm is to explain 
multiple-balancings in terms of what are perceived to be "true" reactions, namely 
those admitting unique balancings. In fact, any reaction can be written as a ra- 
tional combination of reactions with the reactants (resp. products) being subsets 



of the original reactants (resp. products) though as demonstrated in Example 3.4 
we cannot require the subsets always to be proper. This follows from a simple 
observation: let convQ denote convex hull. 

3.6. Lemma. Suppose conv(X) = conv(P) n conv(P) for finite sets P, P in some 
fixed Euclidean space. Then, there exist subsets Ri C R and Pi C P such that 
the intersections conv(Pi) n conv(Pi) = {xi} are singletons and any point x € 
conv(X) is a convex combination of the points X{. We cannot require Ri,P% to be 
simultaneously proper. 

Proof. We may assume that X is the set of generating vertices for the intersection 
conv(i?) Plconv(P). In general, for each Xi € X, let Ri,Pi be minimal subsets of 
R, P respectively such that xi € conv(Pi) n conv(Pj). If Xi is an interior point of 
either conv(P), conv(P) then a set of generating vertices for one of the two hulls 
may be needed. Hence, we cannot require Pi, Pi to be necessarily proper subsets. 
If Xi is a boundary point for both conv(P), conv(P), then it is the intersection of a 
line segment of one hull and a face of the other hull. In either case the intersections 
conv ( Ri )flconv( Pi) are singletons. Since the X4 generate conv(JT), we are done. □ 

Taking P, P to be the sets of reactants and products, respectively and conv(X) 
to be the intersection polytope for the reaction, the we obtain the following: 

3.7. Proposition. Any chemical equation can be written as a rational combination 
of chemical equations admitting unique balance, with the reactants (resp. products) 
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being subsets of the original reactants ( resp. products ) though we cannot require the 
subsets always to be proper. 

Proof. Refer to the proof of the preceding lemma. All that remains to show is that 
each Xi encodes an unique balance. Note that any point of a finitely generated 
convex set is in the interior of a simplex generated by some (necessarily) minimal 
subset of the generating vertices. This simplex need not be a face of the polyhedron. 
If Xi is an interior point, then the preceding observation shows that Xj is in the 
interior of some simplex. If boundary point, then it is the intersection of 

a line segment and a face so that again by the opening observation, we may think 
of Xi as the intersection of a line segment with a simplex. In either case, since the 



vertices of a simplex are linearly independent, by Equation (2.18), Xi represents an 



unique balance. □ 

This thus justifies the chemical paradigm of expressing chemical equations as 
mixtures of those admitting unique balance. 

4. Balancing Equations Involving Charged Species 

If charged species are involved, there are several ways we may procede: we 
may introduce (a) an additional component to index the charge or (b) fictitious 
"spectator ions" that allow us to balance the equation as if only neutral species 
are present. One important class of chemical reactions involving charged ions are 
oxidation-reduction reactions. Though these may be balanced by either of the 
two preceding methods, there is a commonly used "half-reaction method" that 
merits separate mathematical investigation since it is structurally different from 
the preceding methods. 

4.1. Component Indexing Charge. We index charges as an additional compo- 
nent so that if n elements are involved, then the species vectors are in M. n+1 now. 
We use the convention that the charge is always indexed last and note that it 
may be negative. For example, the species Ca 2+ would have species vector (1,2) T 
whereas Cl~ would have (1,— 1) T . However, balancings still correspond to equali- 
ties of non-negative rational linear combinations, and almost everything above for 
both the algebraic and geometric methods remain valid. Evidently, the codomain 



of R,P will be so that the commutative diagram (2.17) should be altered 

accordingly, but the subtlety is that since we now allow negative entries in the last 
component, we no longer have a canonical choice for the slicing hyperplane. 

If both lone +/— charges (i.e. not bound with an element) are present, then no 
hyperplane always works because any hyperplane intersecting both the positive and 
negative a;„+i-rays contains the a;„+i-axis and hence passes through the origin. Of 
course, one can assume that the lone charges have been "canceled" so that only one 
kind remains. Alternatively, since in chemistry lone positive charges are denoted 
H + , we may assume only lone negative charges are possibly present. 

Recall, that the choice yV., Xi = 1 is canonical if charges do not appear by 
themselves without elements, e.g. free electrons e~ . When only one kind of free 
charge is involved (either positive or negative), we may be tempted to consider 
as a canonical hyperplane: ±x„+i + Xa=i x i = 1> where the sign depends on the 
parity of the free charge. Unfortunately, as the following example illustrates, this 
hyperplane will not necessarily intersect both the "reactant" and "product" cones: 
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4.1. Example. Consider the reaction: 

H N0 2 -> N0 2 +H+ + e~ 

Because a lone negative charge is a specie, we use the hyperplane — £4+X)i=i x i = 1 
(remember that charge is indexed last). However, the ray corresponding to H + is 
of the form (fc, 0, 0, k) for real k > and thus does not intersect the hyperplane. 

Even without a canonical slicing hyperplane, there still exist (many!) hyper- 
planes which slice both "reactant" and "product" cones of a given reaction. If we 
insist on working with polytopes instead of cones, we need only slice with any such 
hyperplane, albeit non-canonical. 

4.2. Spectator Ions. One method chemists utilize for balancing equations involv- 
ing charged species is to introduce fictitious spectator ions (e.g. Q + ,X~) which 
are formal symbols bound with the charged species so that the equation no longer 
involves charged species. After a balance is found, the spectator ions are deleted. 
The following example illustrates this approach: 

4.2. Example. Consider the unbalanced equation: 

Zn + NOz + H+ -> Zn 2+ + NH+ + H 2 (unbalanced) 

Using spectator ions Q + ,X~ we instead balance 

Zn + N0 3 Q + HX ZnX 2 + NH A X + H 2 + QX 

Note the additional inclusion of QX among the products so that Q will be present 
on both sides. This yields an unique balance: 

AZn + N0 3 Q + IQHX -> AZnX 2 + NH 4 X + 3H 2 + QX 

Hence, the desired balance to our original equation is 

AZn + ATO" + 1QH + -> AZn 2+ + NH+ + 3H 2 

The geometric explanation is simple: Let £ be the non-negative span of the 
elements involved in the reaction. Then the image of the species simplex lies in 
£ x R where the last coordinate indexes charge. Thinking of K as the union of the 
non-negative and non-positive rays, we have £ xl = £ X (K>oUK<o) as sets (simply 
by twisting one of the rays orthogonally) . Identifying the latter two components as 
the non-negative span of X, Q respectively, we see that the codomain may be taken 
as the simplex spanned by the elements together with the spectator ions. Each 
specie is represented by a point in either the £Q, £X-hyperplanes or the point QX 
in the QX-plane, though the reactant and product cones are not restricted to these 
hyperplanes. In fact, the introduction of QX is necessary in situations where the 
reactant or product polytope lie in different boundary faces of the orthorant of 
£ x M>o x M<o or only one does while the other is in the interior. See Figure [8j 

4.3. Half- Reaction Method. One important class of chemical reactions are called 
oxidation-reduction reactions. They are characterized by changes in the oxidation 
states of elements and often involve H + ,OH~ , H 2 as species. In fact, these 
species are sometimes not explicitly specified until an equation is balanced. For 
chemical details, see any general chemistry textbook (e.g. Qlmsted|19j). One com- 
mon approach to balancing oxidation-reduction reactions is the "half-reaction" or 
"ion-electron" method. This method first divides the reactant and product species 
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Figure 8. Geometric Explanation for Spectator Ions 
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Figure 9. Geometry of Balancing Redox Half- Reaction: Steps (3) 
and (4). 



into two "half-reactions" corresponding to an "oxidation" and a "reduction" reac- 
tion. This selection is determined by chemical considerations. The half-reactions 
are separately balanced using H + , OH~ , H2O, e~ and then linearly combined so as 
to cancel unbound elections e~ . The chemical considerations assure that unbound 
electrons in the two half-reactions will occur on opposite sides of the equation and 
hence are amenable to an unique cancellation via linear combination. Thus, the 
result of the half-reaction method is unique iff each half-reaction admits an unique 
balance. 

Each half-reaction is balanced using the following recipe: 

(1) balance all elements except H, O, and the charge 

(2) balance oxygen O with H 2 

(3) balance hydrogen H with H + ions 
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(4) (for basic solutions): cancel H + using the equality H 2 = H + + OH 

(5) balance charge with electrons 

The geometry of steps (3) and (4) are illustrated in Figure [9] Symbols X, Y repre- 
sent the complexes (i.e. each individual side of the equation) after step (2), shown 
by their projections in directions 0,H only. The diagram shows the case where 
complex X is deficient in H and hence balanced with H + resulting in V. Symbol 
W denotes water H2O and Z is the result of (4) (recall both charge directions are 
projected out). Line segments H + W and VZ are parallel. The remaining step 
(5) balances the two lifts of Z (as reactant and product complexes) from the face 
spanned by the elements into its join with the charge (i.e. spectator ion) directions. 

Notice that the above recipe provides an unique balance to a half-reaction iff 
the projection onto the orthogonal complement of the subspace indexing H, O, and 
charge admits an unique balancing. All steps subsequent to the first are uniquely 
determined and provide a sequence of lifts of the partial solution into the subspaces 
indexed by O, H , and charge, respectively. 

The following example illustrates the "half-reaction method" : 

4.3. Example. [19, Ex 17.7, p827] Consider the following reaction in basic solution: 

Au + 2 + CN~ -> [Au(CN) 2 ]~ + H 2 2 (redox reaction) 
By chemical principles one identifies that gold Au is oxidized by cyanide ions CN~ 

(4.4) Au + CN~ -> [Au(CN) 2 ]~ 
and oxygen gas 2 is reduced to hydrogen peroxide H2O2: 

(4.5) 2 -> H2O2 



The oxidation half-reaction (Equation 4.4) is readily balanced as 



Au + 2CN~ -> [Au(CN) 2 }- + e 



and the reduction half-reaction (Equation 4.5 ) as 

2e" + 2H+ + 2 ^ H 2 2 2e~ + 2H 2 + 2 -> H 2 2 + 20H~ 

since the reaction is taking place in basic solution. Hence, the overall reaction is 

2Au + 4CN- + 2H 2 + 2 -> 2[Au(CN) 2 r + H 2 2 + 20£T 

However, the directness of the recipe belies a subtlety: the intersection cone 
associated to the unbalanced reduction half-reaction 

H 2 + 2 + e~ -> H2O2 + 20H- 

has dimension 2, so that both this half-reaction and the total reaction admit multi- 
ple balancings. The subtlety lies in the fact that for basic solutions the recipe does 
not balance using H2O, OH~ directly as in the preceding equation, but instead uses 
H + and then converts via H2O = H + + OH~ . For the preceding example, this 
choice constrains the balanced reduction reaction to have H 2 0, OH~ in a 1 : 1 ra- 
tio. In effect, it slices the intersection cone for the half-reaction with the hyperplane 
defined by equal coefficients for H 2 0,OH~ . Indeed, the one-parameter family of 
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balancings for the reduction action yields a one-parameter family of balancings for 
overall reaction: 

(4.6) (1 - t)Au + 2(1 - t)CN~ + f^P) H 2 + 0^) 2 ->■ 

(1 - t)[Au(CN) 2 }- + tH 2 2 + (1 - t)OH~ 

for i e [0, 1]. The output of the "half-reaction" method corresponds to t = 1/3. In 
particular, the reaction: 

(4.7) 2Au + ACN~ + 4H 2 + 20 2 -> 2[Au{CN) 2 }- + 3H 2 2 + 20R- 
from using the reduction half- reaction (t — 3/5): 

2e- +AH 2 + 20 2 -> 3H 2 2 + 20H~ 

is not obtainable via our two chosen half-reactions since H 2 0,OH~ are not in a 
1 : 1 ratio. 

One may wonder if Equation |4. 7| were obtainable using the "half- reaction" recipe 
relative a different set of half-reactions. This is not the case: for our redox reaction, 
we may classify half-reactions by whether [Au(CN) 2 ]~ is a product. If so, then 
Au, CN~ need be present as reactants to admit a balance. The only other possible 
product combination is H 2 2 alone; since Au, CN~ cannot be included among the 
reactants without requiring [Au(CN) 2 ]~ among the products, we have only four 
possible half-reactions: 



Unbalanced 


Balanced 


Au + CN- -> [Au{CN) 2 ]- 




Au - 


f 2CN- -> [Au(CN) 2 ]- + e~ 


Au + CN~ +0 2 ^ [Au{CN) 2 ]- 


3e~ 


+ 2H 2 - 


VAu + 2CN- -> [Au{CN) 2 ]- + AOH~ 


Au + CN- -> [Au{CN) 2 ]- + H 2 2 


Au - 


- 2CN- - 


f- 20ff- -> [Au(CiV) 2 ]- + # 2 2 + 3e- 


2 -> H 2 2 




2e~ 4 


- 2F 2 + 2 -> i/ 2 2 + 20ff" 



where each has been balanced using the recipe for balancing half-reactions. One 
checks that no pair of half-reactions yields Equation |4.7| Hence, the half-reaction 
method cannot necessarily find all possible balances. 

Recall that the half-reaction method outputs a balancing iff there exist two half- 
reactions for which step (1) above is realizable, i.e. the projection onto the orthogo- 
nal complement of the subspace indexing H, O, and charge admits a balancing. This 
orthogonal complement is equivalently the span of all non H, O-elements. Suppose 
a balancing for the overall reaction exists; then the projection onto the aforesaid 
complement is non-empty. Now, denote the set of "reactant" and "product" species 
by R,P respectively. Hence, step (1) is realizable iff there exists decompositions 
R = Ri U R 2 and P = P x U P 2 such that conv(i?j) n conv(P 4 ) ^ for i = 1,2. 
Here, conv() denotes the convex hull of sets. The following example shows that we 
cannot require the subsets to be disjoint: 

4.8. Example. |14j Consider the following unbalanced oxidation-reduction reaction: 

P 2 h + Pa + H2O -> PH 4 I + H 3 P0 4 (unbalanced) 

If we order P before I in indexing elements, under the projection onto the span 
of P,J the coordinates for the reactants are (2,4) T , (4, 0) T and that for products 



(1, 1) T , (1,0) . A glance at the polytope diagram (see Figure 10 1 shows that we 
cannot require both i?i,i? 2 and P\,P 2 to be disjoint pairs. 
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P«,H 3 P0 4 PH 4 I 




Figure 10. Cannot Partition Reactants and Products into Dis- 
joint Pairs 

The following lemma assures us that otherwise existence follows readily: 

4.9. Lemma. Let conv(X) — conv(i?) n conv(P) where conv(X) has dimension 
c < d for non-empty finite sets R,P,X C M. d where we may assume each set 
individually to be convex linearly independent. 

Let x £ conv(X). Then, there exist decompositions R — R1UR2 and P = P1UP2 
such that defining Hi — conv(_R;)nconv(p) fori — 1,2 we have x € conv(iJiU-ff2)- 
Furthermore, if c is positive then we may require Hi to be non-empty. 

Proof. We first provide an intuitive argument. Since conv(X) is convex, excepting 
the trivial case when conv(X) is a singleton, we may see x by peering through any 
boundary facet of conv(X). By shifting ourselves slightly, we can place x in the 
interior of a background facet of conv(X). Then, x is in the convex hull of these 
two facets of conv(X). 

To formalize the above argument, we place conv(X) in a c-ball (where c is the 
dimension of conv(X)) such that the center of the ball is in the interior of conv(X). 
Taking the radial projection outwards from the center point yields a cellular decom- 
position of the boundary sphere. Choose any cell. By appropriate choice of point 
p within this cell, the line through p, x meets the sphere at another point q interior 
to some other cell on the sphere (here we use the fact that conv(X) is positive- 
dimensional). These two cells identify two facets H\,Hi of conv(X) such that their 
convex hull includes x. Since each face of conv(X) is the intersection of some faces 
of conv(i?), conv(P) there exist subsets Ri, Pi such that Hi = conv(Ri) nconv(Pj). 
Finally, note that since Hi are facets of conv(X), at least one of Ri, Pi is a proper- 
subset for each i = 1,2. □ 

4.10. Remark. In Example |4.3[ the projection onto the complement of H, O and 
charge yields only the species Au, CN, Au(CN) 2 . The projected reactant and prod- 
uct polytopes meet at a singleton. Hence, we are not assured two non-trivial half- 
reactions in the space spanned by Au, C, N . 

4.11. Proposition. If an oxidation-reduction reaction admits a balance, then there 
exist decompositions of the reactants and products, each into two (not necessarily 
disjoint) groups so that application of the half-reaction method yields a balance 
relative the given reactants and products. By appropriate choice of decompositions, 
this balance can be made to reflect the elemental proportions of all non-O, H elements 
in a given initial balance. 

Proof. As mentioned previously, we need only show that we can decompose the 
reaction into half-reactions for which step (1) of the recipe is realizable. Since the 
reaction admits a balance, the intersection cone remains non-empty after projecting 
out the H, O and charge-directions. Hence, the projections of the reactant and 
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product polytopes satisfy the hypotheses of the lemma. We let Ri , Pi be as assured 
by the lemma. Then, Hi = Ri n Pj is non-empty so that step (1) is realizable for 
each half-reaction. □ 

The half-reaction method has been criticized by some chemists because some- 
times the choice of oxidation/reduction reactions is unclear and even chemically 
fictitious [15]. As discussed above, the method has mathematical drawbacks in that 
the recipe may obscure the presence of multiple balancings and may not capture 
all possible balances. However, at least the preceding proposition assures us that if 
a balancing exists, then the "half-reaction method" can generate a valid balance, 
though the half-reactions may be chemically fictitious. 

5. Reaction Mechanisms and Stoichiometric Restrictions 

In this section we discuss the relationship between reaction mechanisms and 
stoichiometric restrictions. Given a chemical reaction, a reaction mechanism is a 
sequence of steps that together elucidate how the reactants of the overall reac- 
tion are converted into the products of the overall reaction. Each step is called 
an elementary chemical reaction and represents the most fundamental molecular 
transformations. 

5.1. Prom Mechanism to Linear Stoichiometric Relations. A given reaction 
mechanism has implications for observed stoichiometrics of the involved species as 
illustrated in the following example: 

5.1. Example. [22] Let Si,S 2 ,S 3 be species and consider the following reaction 
mechanism: 

Si — > S 2 + S3 
Si + S 2 -> 2S 3 

Letting [Si] denote (time-dependent) concentrations, the theory of chemical kinetics 
implies that 



d[S- 



dt 
d[S 2 ] 

dt 
d[S 3 ] 

dt 



— — X\ — X 2 
=X\ — X 2 
=X 1 + 2X 2 



where the Xi depend only on the i-th elementary reaction in the mechanism. Elim- 
inating the Xi yields the single relation: 

* t (3[S 1 ] + [S 2 ]+2[S 3 ]) = 
which constrains the observed stoichiometries. 

We may formalize this as follows: Order the species involved and let A be a 
matrix with columns indexed by the steps of the mechanism and the rows indexed 
by species. We will follow the convention in the literature of negative entries for 
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reactant species and positive entries for product species. For the preceding example, 
we would have 



Then, elimination of the Xj is tantamount to finding elements of NS(N T ) which we 
may think of as the space of stoichiometric restrictions. Each vector in this subspace 
defines a linear relation on the concentration of species that integrates to a stoi- 
chiometric restriction which in the chemical literature is called a mass-conservation 
equation. (The terminology arises from the assumption of constant volume during 
the course of the reaction.) In particular, mass-conservation equations are defined 
with respect to a given set of reactions. Note that though d[Si]/dt = X,-, since we 
do not assume a particular form of Xi we need not assume mass-action kinetics. 
This shows that for given initial concentrations [Si] , the resulting flow is constrained 
to lie in an affine space linear isomorphic to M = CS(iV), the orthogonal comple- 
ment to NS(N T ). Further information can be obtained if we assume mass-action 
kinetics, whence we obtain an ODE system of polynomials. The literature on these 
systems is vast (see for example either Horn and Jackson[13] or Erdi and T6th|9 ]). 

In addition, chemists also define element conservation equations. These are de- 
fined for any given single reaction. Given the elemental proportion matrix M for 
a reaction, since the rows index elements, we may require their conservation in the 
course of a reaction. Thus, the set of all element-conservation equations corresponds 
to the row space RS(M). 

Up to now our discussion has not assumed any elemental composition for the 
species nor need there be one. The existence of such implies a (non-negative) el- 
emental proportion matrix M such that MN = because each reaction indexed 
by N is balanced. In particular, we may assume the rows of M index only ele- 
ments present so that the row space RS(M) must contain a positive vector. Since 
NS(iV T ) D RS(M), a necessary condition for the existence of M is that NS(7V T ) 
contains a positive vector. In the chemical literature, such matrices N are known 
as conservative |18j . The sufficiency of this condition is readily checked. Further- 
more, given an overall reaction and a reaction mechanism for it, since the rows 
of M index elements, the relation RS(Af) C NS(7V T ) shows that any elemental- 
conservation equation is a mass-conservation equation but not vice-versa (see also 
Smith and MissenrTS]). Letting X,S,£ denote the spaces spanned by elementary 
reactions, species, and elements respectively and thinking of N, M as linear maps, 
we have 



and the discrepancy between mass- and element-conservation equations is described 
by the homology (ker M) / (im N) at S. 

In a reaction, there are usually short-lived species called intermediates that are 
difficult to observe. Under a steady-state hypothesis, or if we consider only com- 
pleted reactions, or just assuming that they cannot be detected, the space of ob- 
served linear dependencies is instead the orthogonal projection of NS(iV T ) onto the 
span of observable species. We denote this space ir^. NS(iV T ). 




X -> S -> £ 



5.2. Further Conclusions from a Given Mechanism. In this subsection we 
show how several chemically interesting problems lead to the problem of counting 
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lattice points in polytopes, a much-studied and active research area in combina- 
torics. 

5.2.1. Reactions Consistent with Mechanism. Let N be an integer matrix associated 
with a collection of (elementary) reactions, and write AT = CS>o(-/V) for the cone 
determined by its non-negative span. Let K be the set of known (or observable) 
species and IC their span. For a chemical equation c 6 IC we will say that c is 
consistent with N or that N is a mechanism for c provided c £ Af . 

For a given N, the set of consistent chemical equations is thus the set of lattice 
points in AfnJC modulo dilation, or alternatively the points of the projective subset 
Pq(jV n IC). Suppose the cone Af (1 JC is non-trivial. Then, Pq(Af n JC) is either 
a singleton representing an integer lattice point or contains infinitely-many such 
representatives. One can justify the latter claim by noting the infinitude of rational 
vectors in any non-trivial projective neighborhood of an integral vector. 

Alternatively, via Ehrhart theory the number of lattice points in integral dilations 
of a rational ei-polytope P are counted by a quasi-polynomial Ep called the Ehrhart 
quasi-polynomial (see [7J, [5]). That is, there exists some N > and degree d- 
polynomials /o, • • • , /jv-i sucn that Ep(t) — fi(t) for t = i mod N. Lattice points 
in the relative interior are counted by another quasi-polynomial Ep[t) where the two 
quasi-polynomials are related by Ep-(t) = (—l) d Ep(—t) where d = dim P. Hence, if 
d > 1, then E P grows at least quadratically so that Fq^AFCiJC) contains infinitely- 
many points. Chemically, this implies that there is either an unique chemical 
equation consistent with a given mechanism or there are infinitely many. However, 
since we do not expect to observe overall reactions built from many instances of 
each elementary step, in practice we see but finitely-many consistent equations. 

Suppose we want to consider overall reactions containing at most t species (count- 
ing multiplicities) . Then, we can count the lattice points within the polytopc formed 
by bounding the cone Af tlK. with the cross-polytopes Cf = {z e R d \ \ z i\ = £}■ 
This is computationally feasible since there exist polynomial-time algorithms to 
compute the number of lattice points in polytopes. See deLoera[7j or Barvinok[2],[3] 
and the programs LattE[6] and BARVINOk[24]. (The output of the latter is techni- 
cally a piece-wise step-polynomial, a generalization of quasi-polynomials.) However, 
we should keep in mind that the counts from lattice point enumeration in S may 
include multiples of a given balance. 

In light of the above discussion, we may be tempted to define a partial order 
on the interior points of the cone Af fl K via the total number of known species 
involved (counting multiplicity), namely YIk Sj. However, there exist lattice points 
on arbitrarily large cross-polytopes such that their convex hull with the origin does 
not include any interior lattice point. For example, taking i,j, k unique indices, let 
x = tei — ej and y = tej — for positive integer t. The convex hull conv({0, x, y}) 
contains no interior point. Thus, relative the ^-norm there may not be a minimum 
consistent reaction. 

5.2. Remark. We may also want to consider consistent reactions with mechanisms 
involving no more than a fixed number of reactions. We then would have to work 
in reaction space X by counting lattice points in the intersection of cross-polytopes 
and A>o n N^ 1 (Af n IC) where X>q is the non-negative orthorant in X . 

5.3. Remark. This observation on lattice points of cross-polytopes also applies to 
studying the polyhedron Q C S that parametrizes balances of a given chemical 
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equation. Hence, using the sum of the coefficients in a balance may not be sufficient 
to single out a unique "simplest" balance. However, we may use the Ehrhart theory 
to estimate the minimal total number of species involved in the overall reaction 
(counting multiplicity) by intersecting Q with cross-polytopes C t s in S . 

5.4. Remark. We express Craciun and Pantea's result on confoundable reaction 
networks using our framework. Two reaction networks are confoundable if they 
yield the same kinetic differential equations under the assumption of mass-action 
kinetics. 

Let N be the equations of a reaction network expressed as a matrix (columns 
indexing reactions, rows species), S the vector of specie concentrations, and fe, the 
forward rate-constant for the i-th reaction. Let x be a vector of dimension equal to 
the number of elementary reactions in the mechanism, and suppose that the entries 
of x depend on a function of the reactant (i.e. source) complexes alone. Then, 
dS/dt — iVdiag(fcj)x. We may assume the components of x and the columns of 
iV are sorted by reactant complex so that iVdiag(fcj) has only as many columns as 
distinct source complexes and the dimension of x is the number of distinct reactant 
complexes. (Just add together columns of N diag(fci) indexed by reactions having 
identical reactant complex so each column of N diag(fc^) is a cone on N parametrized 
by the fcj of reactions with that reactant complex.) Thus, a mechanism associates to 
each reactant complex a cone. Then, two networks yield the same kinetic differntial 
equations provided they have identical set of reactant complexes and for a fixed 
reactant complex, the cone so indexed for each of the two mechanisms intersect 
non-trivially. This is Craciun and Pantea's result, though without assuming mass- 
action kinetics. 

5.5. Example (Decomposition of Azomethane) . We consider the decomposition of 
azomethane: 

(5.6) 5C 2 H 6 N 2 -> 3N 2 + CH 4 + C 2 H 6 + C 3 H S N 2 + C 4 H 12 N 2 

A reaction mechansim is given by Missen and Smith [18 which involves intermedi- 
ates that we will label X, Y, Z. Under the ordering 

C 2 H 6 N 2l N 2 ,CHi, C 2 H 6 , C 3 H 8 N 2 ,CiH 12 N 2 , X, Y, Z 

the mechanism has associated matrix 
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To count consistent overall reactions where no more than a fixed number of 
species are present in the equation, we consider the intersections (cone(0, N) H/C) n 
C\ where Cf are cross-polytopes. For example, the intersection with C| is a 6- 
polytope in S = M 9 having 14 vertices. By barvinok we find that this polytope 
has 35 lattice points so that there are 35 consistent overall reactions where there 
are no more than 6 species (counting multiplicity). 
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3-f 

NS(M)» • RS(M) 



NS(N T ) 

Figure 1 1 . Relations among subspaces of S 

5.2.2. Algebraic Representations of a Reaction. Let c £ K, be a given chemical 
equation with mechanism N . We will call any non-negative integral solution to 
Ax = c an algebraic representation for c. Algebraic representations do not nec- 
essarily correspond to chemical pathways within a given mechanism because the 
presence of intermediates in elementary reactions means that the order of reactions 
in a mechanism is important. For example let I±, I2 be intermediates, and suppose 
the following pair of equations: 

S 1 + I 1 = S 2 + h s 3 + h = S A + h 

comes from some mechanism. Then their sum is an algebraic representation for 
S1 + S3 = S2 + S4, but the two reactions alone cannot constitute a chemically viable 
mechanism unless one of the intermediates is present at the start of the reaction. 
Further discussion of this restriction is postponed until Subsection |5.4| For now we 
note that the number of algebraic representations is given by the number of non- 
negative lattice points in the polyhedron determined by Ax = c. As mentioned 
above, this is computationally feasible to determine. This number is finite iff NS(A) 
has no non- negative element, or equivalently that is not in the convex hull of the 
columns of A. If infinite, we may still study representations involving no more than 
i-steps in total by slicing the non-negative solution cone of Ax = c in X with the 
hyperplane yV, Xi =t where x = dim X the number of columns of A. 

The case of unique algebraic representation corresponds to polytopes {x|Ax = 
c} with unique interior lattice point. Because we are working over Z, this may 
happen even if ker JV^O. For example, the polytope corresponding to the intersec- 
tion of the hyperplane {{x,y)\x + y — 2} with the non-negative quadrant is a line 
segment with an unique interior lattice point. Polytopes in K 2 with unique interior 
lattice point have been classified up to unimodular equivalence in Rabinowitz |20j . 
but seems unknown in higher dimensions. Already, there exist lattice simplices in 
R d with unique interior lattice point but having more than 2 2 boundary lattice 
points (see Zaks, Perles, and Wills [2*5]). 

5.8. Example. Again, consider the mechanism for the decomposition of azomcthanc 
given in Example |5.5| The convex hull of the columns of A determines a convex 
5-polytope in S = E 9 . Since this polytope does not contain 0, there are only finitely- 
many non-negative solutions to Ax = c for any c 6 S. Using c = (—5, 3, 1, 1, 1, 1) T 



from the reaction (5.6 1, we find that the polytope parametrizing non-negative so- 
lutions is actually just a point, namely (3, 1, 1, 1, 1, 1) T . Thus, this is the unique 
non-negative integral combination of elementary reactions of the mechanism N 



yielding the overall reaction (5.6) 
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FIGURE 12. Relation among subspaces of K. The possibility of 
non-trivial intersection between O and Z — (tt/c RS(M))- l has been 
suppressed in the Diagram. 



5.3. Mechanisms Consistent with Observed Relations. We now consider an 
inverse problem. Suppose we know all the species (including intermediates) involved 
in a reaction with unknown mechanism; denote their span by S. A set of linear 
dependencies are observed which span a space O (obtained via multivariate linear 
regression of data, for example) . What can we say about mechanisms -/V compatible 
with these restrictions and in particular of the image space Af = imiV? (Note here 
Af represents a linear space contra the cone in the preceding subsection.) We know 
that for a given chemical equation c, the moduli of mechanisms for c is given by 
the set of all integer matrices N such that c £ Af. 

Figures [TT] and [T2| summarize the relations among the subspace now to be defined. 
Let K,U be the sets of known and unknown (intermediates) species, respectively 
and S = K U U. Write JC,U,S for the spaces they span and equip S with the 
canonical inner product. Then S = K, _L Li. Now, write kerM = Af JL H where we 
think of H as the homology at A S 4. Then, S = Af ± U -L RS(M) so that 
H ± RS(M) = NS(7V T ) = NS(7V T ) _L tt w NS(iV T ) 

where tt/c^u denote orthogonal projection onto ]C,U respectively. Write O for 
the space of observed dependencies. We will not assume that this constitutes all 
possible linear dependencies, hence we have 

(5.9) ix K U _L tt k RS(M) = ttk; NS(7V t ) = O _L W 
for some subspace W . In particular, we see that 

(5.10) dim ttk; NS(7V t ) > dim(C + tt^ RS(M)) 



where the right-hand-side may be computed without knowledge of any intermedi- 
ates. If we specify intermediates, we can do much better. From Equation (5.9| we 
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have 

(5.H) [O + W] C [ttjcH] C ^ 

7Tk; RS(Af) 

and since K. = tt/c NS(Af) + 7r^ RS(M) we have 

< 5 ' 12 > ».NS(M) N n S( "'R S (M) " i^(M) " (* RS ( M » 1 - Z C « NS < M > 
so that we must choose a subspace X such that 

(5.13) proj z 0CX=: P roj z H 

where we choose a proper superset when we believe that W ^ 0, i.e. that do not 
represent all possibly observable stoichiometric restrictions. Regardless, we have 

(5.14) 7r K -H = X + ( 7 r ;c NS(Af)n7r yc RS(M)) 

D proj z + (Tryc NS(M) n tt k RS(M)) 

Since H C NS(Af), to lift 7rycH back up to NS(M) we need to study the kernel 
of the restriction of n/c to NS(M). We denote this map ip — tt/c : NS(M) — » /C. 
Writing M = (Mr- Afy) where Mjc,Mu are the submatrices of M indexed by 
known and intermediate species, respectively, we see that NS(M ) is the centralizer 
of the maps Mijc^tc^ —Miu^u so that ker</? = in NS(A/[/). Thus, 

(5.15) H = tt^X + ^{-KK NS(M) n tt k RS(M)) + i M NS(M t /) 

D tt" 1 proj z + 7r -1 (7rx NS(M) n ttk RS(M)) + i u NS(M £/ ) 

In the case NS(iV T ) = 0, we have X — proj z so that the possible choices for H 
(and hence AT since AT = (H © RS(M)) J -) are parametrized by subspaces of 

(jtjc NS(Af) n 7r c RS(AO) -L iu NS(A%) 
with dimension bounded by dim(0/Z). We illustrate these ideas with two examples: 

5.16. Example. Again, consider the decomposition of azomethane given by Equation 
(5.6) of Example 5.5 above. With the species ordered as in that example and 
ordering the elements C, H, N in that order, the associated elemental proportion 
matrix is 

/2 1 2 3 4 1 2 3\ 
M = 60468 12 359 
\2 2 2 2 2 2 y 

where the intermediates X, Y, Z have molecular formulae CH3, C2H5N2, C3HQN2 
respectively. Suppose the space c fC = R 6 of observed stoichiometric restrictions 
has a basis 

(-1,-3, 1,2, 0,1) T (2,4,0,-2,0,0) T (0,0,-l,0,l,0) T 

Let Z = (ttjcRS(M)) ±k . We compute that dim(0/Z) = 0. Thus, assuming 
includes all possible stoichiometric restrictions, we may set H = 0. 

Also, we find that dimT: K: RS{M)nn K: NS(M) = 3 and that NS(A%) = so that 
the possible Af are parametrized by 

Af = (span(s) + RS(M)) 1 - s e tt- 1 (ttkRS(M) n tt k NS{M)) 
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5.17. Example (Oxidation of Formaldehyde). Missen and Smith provide a mech- 
anism in |18j for the oxidation of formaldehyde. If we order the elements by 
H,C,0, Co, + and the species by 

CH 2 0, Co 3+ ,Co 2+ ,H+, CO, H 2 0, CH 2 2 , H 2 ,CH 2 OCo 3+ , CHO, OH, H 

where the last four species are intermediates, we have the associated elemental 
proportion matrix 



(5.18) 
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has associated matrix 
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so that the space O = tt k NS(A T ) c K, = R 8 has a basis 

(1,0, 0,0, 1,0,0, if (0,0, 0,1, -2, 0,0, Of 
(0,1, 0,0, 2, 0,0, Of (0,0, 0,0, 0,1,0, If 
(0, 0, 1, 0, -2, 0, 0, Of (0.0, 0, 0, 0, 0, 1, -if 

Suppose we started knowing only the species and the above space O. Set Z = 
(tt/c RS(M)f' c and let proj z denote the orthogonal projection onto Z. We find 
that dim(0/Z) = 1 where 



proj 2 = span I 1 



-12 12 12 



23 23 23 



-49 

~23~' 



-26 26 -29 

~23~' 23' ~23~ 



Assuming all stoichiometric restrictions are included in O, i.e. O = 7TjcNS(^V t ), 
we have tt/cH — proj z O. A basis for tt/c NS(Af) n tt/c RS(M) is given by 



(2,-1,0, 0,0,0,0,2f 
(-1,0, 0,0, -1,1,0, of 



(1,0, 0,0, 1,0,1, of 
(0,2, 1,1, 0,0,0, Of 



and this represents the uncertainty in determining ttjcH- One checks that ttjc is 
injective on NS(M) so that proj z O + (7r yc NS(Af) H n K RS(M)) lifts uniquely to 
NS(M). Hence, compatible subspaces N = (H _L RS(M))- 1 - are parametrized by 
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the 4-dimensional space ttk. NS(M) n tt/c RS(M). Setting tt^H — proj z O yields 

T 



H 



-14 26 -26 -49 12 -12 



1 



-29 

~23~' ~23~' 23' ~23~' ~23~' 23' ~23~ 
However, using the mechanism N given above by Missen and Smith yields 



40 26 
Spanl 23'23' 1 ' 1 ' 



H = span 1, — — -, — -, — , 



101 17 -43 



219 315 V 



-75 -5 7 

64 '32' 32' 3*T'128'4' 128 '64' 64 ' 128 '128,/ 
The difference of their projections into K. is 



761 



-1573 -9 -9 -4379 13249 4435 1285 Y 



12928' 6464 '404' 404' 12928 ' 12928 ' 6464 ' 3232/ 
which one verifies is in the subspace ttjc NS(M) n tt/c RS(M). 



The dimensions of important subspaces in Examples |5.16| and 5.17 
rized in the table below: 



are summa- 



Space 


Ex 5.17 


Ex 5.16 


NS(M) 
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RS(M) 
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7T/C NS(M) 
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TTyc RS(M) 
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tt k NS(M) Httk; RS(M) 
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e> 
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G/^RS(M) 


1 





tt w RS(Af) 


5 


3 



5.4. Prom Subspaces JV to Mechanisms. Two tasks remain in obtaining reac- 
tion mechanisms consistent with observed stoichiometric restrictions. The first is 
convex-geometric: we seek spanning sets in Af whose positive hull contains the over- 
all reaction. This ensures algebraic representability. The second is order-theoretic: 
these algebraic representations must then be subject to the constraint on the ap- 
pearance of intermediates discussed in Subsection |5. 2. 2| 

The first task seems computationally expensive: Elementary reactions have usu- 
ally no more than two (or at most three) molecules of a given specie so that we can 
restrict attention to those lattice points in AfP\ [—3, 3] s where s = dim S having both 
positive and negative entries. However, this grows exponentially with the number 
of species involved. Perhaps some savings may result from working projectively 
and this may also be helpful enumerating spanning sets, but on this the author is 



uncertain. Nevertheless, as explained in Subsection 5.2.2 algebraic representatives 
may be obtained from the lattice points. 

The second task is more tractable. Let X C X be a finite set of reactions, 
u the number of distinct intermediates among the reactionsin X, and A" the u- 
dimensional simplex (as an abstract simplicial complex). We may assume there 
are no redundacies in X. We label the vertices K,U\, . . . ,U U . Define set maps 
R,P : X —> A" mapping a reaction to the simplex indexed by its reactants and 
products respectively, with knowns collapsed to K. For example, if we have knowns 
Ki 1 K 2 and intermediates Ui,U2, the maps R,P take the reaction K i + Ux — > 
K 2 + U 2 to the 1-simplices [K, Ui] and [K, U 2 ] respectively. 
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Figure 13. Precedence relations among the reactions in the mech- 



anism for the decomposition of azomethane. See Example 5.20 and 
subsequent discussion. 



Let V(-) denote the set of vertices of a simplicial complex, and for each set A let 
A (A) denote the simplex (as an abstract simplicial complex) with vertices indexed 
by the elements of A. For any non-empty subcomplex Y C A" define 

<p(Y): = (Ao V)(XUPR- 1 (Y)) 

The coface to the vertex K in the simplex ip(Y) indexes all intermediates involved 
(either as reactants or products) in the reactions indexed by X. We set </s({0}) = 
{if}. Since we assume X is finite, the iterates 



0< 



eventually stabilize. In particular the [/-indices of this complex index the inter- 
mediates that can possibly arise through reactions in set X. Note also that if a 
reaction x G X occurs then there is some sequence of reactions x±, . . . , xjy where the 
reactants of xi are known species and the intermediate reactants of Xj are among 
the intermediate products of Xi for i < j. Thus, x G (R^ 1 ° <fi N+1 )({®}) so that 
if a reaction occurs then it is an element of (i? _1 o ^ l )({0}) for some i > 0. The 
converse is readily checked. Finally, note that the sets (i? _1 o v? 4 )({0}) for sucessive 
i > 1 define a partial order on the elements in X the set of reactions. In particular, 
reactions not involved in this partial order cannot occur. 



5.20. Example. Consider again the decomposition of azomethane from Example |5.5| 
where we denoted the intermediates X,Y, Z. Since M = NS(M) we may compute 
that there are 116 non-zero intergal vectors in the intersection NS(M) D [ — 1, l] 9 , or 
equivalently 58 distinct lines. It seems computationally prohibitive to enumerate 
spanning sets for Af among these without further heuristics. 

Nonetheless, choosing as a particular example the set of equations given by the 



mechanism N in Equation (5.7), we have 
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We compute that <p(K) = A({K,X}) and that f 2 (K) = A 3 . The associated poset 
is given in Figure |13| 

Thus, we see that an algebraic representation is order realizable iff it is equal to 
o y l )({0}) for some i. 

5.21. Remark. We have only considered stoichiometric properties of reactions. In 
particular we have not incorporated stereochemical and thermodynamic consider- 
ations and such may help reduce the search space in determining spanning sets for 
computing algebraic representatives. Furthermore, much of the discussion requires 
knowledge of the intermediate species which would need to be guessed by chemical 
intuition. Regardless, it is hoped that the above procedure provides an attempt at 
the computational discovery of reaction mechanisms. 

6. Appendix: Mathematical Justifications for Geometric Viewpoint 

In this section we justify several mathematical assertions from Section [2] First, 
here is the justification of the linear algebraic fact invoked in the discussion relation 
the polyhedron Q and the intersection cone I: 

Let S — CS(^4„ X a) and T = CS(£?„ X b) and consider the n x (a + b) matrix M 
with columns from A, B respectively. Then, S + T = CS(M). 

6.1. Proposition. nullity(M) - din^S" (IT) = nullity(A) + nullity(B). 

Proof. By the Rank-Nullity Theorem for M: 

dim(CS(A/)) + dim(NS(M)) = a + b 

By choosing a basis for S <~) T, separately extending to bases for S, T and then 
invoking inclusion-exclusion: 

dim(5 + T) + dim(5 n T) = dim(S*) + dim(T) 

Subtracting these two equations and using the Rank-Nullity Theorem for A, B 
yields 

nullity(Af) - dim(5 H T) = nullity(A) + nullity(B) 

□ 

Next, we provide geometric justifications for the statements from Section [2] that 
for a given choice of "reactants" and "products" : 

(1) A non-empty intersection polytope necessarily contains a rational point. 

(2) Each rational point in the intersection polytope corresponds to some bal- 
ancing of the chemical equation. 



Algebraic arguments essentially follow from "chasing" Diagram 2.17 
The former follows readily from the following two results: 

6.2. Proposition. The intersection of two polytopes with rational vertices is a 
polytope with rational vertices. 

Proof. Since the polytopes have rational vertices, each supporting hyperplane may 
be chosen to be of the form ^2 niXi — h for rational n^, h. Now, each vertex in the 
intersection polytope is given as the unique intersection locus of certain supporting 
hyperplanes, hence the unique solution of a system TVx = h where N, h have 
rational entries. Thus, each vertex is rational. □ 
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6.3. Proposition. A positive-dimensional polytope with rational vertices contains 
a rational point in its relative interior. 

Proof. Let d be the dimension of the polytope and {x^} any d+1 generating vertices 
whose afHne hull is <i-dimensional. Any convex linear combination with 
U rational and < U < 1 yields a point with the desired properties. □ 

6.4. Corollary. A non-empty intersection polytope necessarily contains a rational 
point. If the intersection polytope is positive-dimensional, we can find a rational 
point in its relative interior. 

Proof. The "reactant" and "product" polytopes have rational vertices. Hence, the 
intersection polytope has rational vertices, and by the preceding proposition, we 
are done. □ 

6.5. Proposition. Each rational point in the intersection polytope corresponds to 
some balancing of the chemical equation. 

Proof. It suffices to show that a rational point in the convex hull of rational points 
is expressible as a rational linear combination of those points. We proceed by 
induction on the number of generating vertices and may assume without loss of 
generality that the generators are a minimal set for the given hull. The base case of 
a line segment determined by its two endpoints is evident. Suppose the statement 
is true for points in the convex hull of at most d rational points. Let p be a 
rational point in the convex hull of d + 1 rational points. We may assume that 
the convex hull is not generated by any proper subset of the d+1 points. Let q 
be one of the generating points. By using a rational frame based at q, we may 
assume the ambient space is R d . The d generating points besides q determine a 
(d — l)-hyperplane H of the form n • x = b where n, b arc rational and n is normal 
to H. Then, n • p and n • q are rational numbers so that by taking dot product of 
n with the orthogonal decompositions 

(6.6) p = s+pn q = t + <?n 

where s, t are in hyperplane H , we conclude that p, q ^ are rational. Finally, 

(6.7) p = s + ^(q-t) = U--t] +^q 

q \ q J q 

a rational linear combination as desired. □ 
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